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1. Introduction 

Entanglement, is a special type of correlation between two or more interacting 
quantum systems [H [2]. Since they are without any classical analogy, generating 
such correlations has become important for the purposes of Quantum Computing and 
Quantum Information [31 H]. Among the experimental demonstrations of entanglement, 
is that with atomic cavity QED systems in the microwave regime [5]. Apart from the 
generation of EPR states [SI [7], a phase gate [HI E] has been realised along with the 
creation of Schrodinger cats [TO] . 

One of the problems encountered when considering atoms traversing a cavity 
resonator, is that of the spatial effects due to the structure of the resonator mode 
As long as the atoms are moving fast [12], then one can safely assume that the 
motion of the atomic center of mass is classical, and utilize the spatial dependence of the 
coupling functions with time dependent pulses [131 [HI |15] . Based on this assumption, in 
recent works we were able to study the adiabatic sequential passage of atoms through a 
cavity [TH [17] . The main feature for the system is the existence of a pure crossing in the 
adiabatic limit, along with a number of possible applications in Quantum Computing, 
including atomic entanglement. 

Here we study the off-resonant limit for the two-atom time-dependent Tavis- 
Cummings Hamiltonian [181 HH], where now the atomic transition frequencies are 
subjected to time-dependent chirps. Although one would expect the adiabatic 
approximation to fail, the main result is that with a proper choice of Stark-shifting 
pulses, the system follows a similar adiabatic evolution as in the resonant limit [16j. In 
addition to this, the use of a chirp is a robust tool for fine tuning the adiabatic phases 
and consequently the output of all possible applications. When taking into account the 
recent developments in the engineering of high quality cavities [201 12T] , we are able to 
demonstrate the generation of a maximally entangled state with very high fidelity. 

The paper is organised as follows: In section [2] the system and the interaction 
Hamiltonian are introduced and after deriving the evolution matrix, we demonstrate 
how a maximally entangled state is formed. Section [3] gives a brief summary of the 
properties of entanglement with respect to the mode structure, and in section jH we 
discuss the importance of spontaneous emission and decoherence. We summarize our 
results in section [51 

2. Time-dependent Tavis-Cummings model 

The system under consideration consists of two atoms sequentially crossing an open 
spherical mirror cavity with velocities f i = f 2 = figure [H The atoms enter the cavity 
with a time delay 2At, where in addition to this they can follow different trajectories 
inside the cavity. The Hamiltonian in the rotating wave picture and within the rotating 
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Figure 1. A pair of atoms with velocities vi — V2 = v enter a Gaussian mode cavity of 
width 2wq with a delay 2Ai. Two Gaussian EM pulses Ej{x) are used to Stark-shift 
the atoms at different times. The y-axis is into the page. 



wave approximation reads {h = 1) 

i=l,2 i=l,2 

where Aj(r) is the detuning of the j— th atom from the cavity mode. 

For a Gaussian mode of width 2wo [TT] , the couphng functions rjj have the following 

form 

^i(^) = ^lexp (-(r + 5)^) , 

(2) 

V2{t) = ^2exp (-(r - 5)2) , 

where in general gi 7^ (72- The dimensionless time r and delay 6 are defined in terms of 
the interaction time cr to be 

t At Wo 

2a 2a V 

In deriving Hamiltonian ([T]) , we are assuming that the center of mass motion is classical, 
and for this to be the case the atoms must be fast to avoid being reflected by the mode 
field m. 



For resonant interactions, i.e. Aj = 0, the adiabatic limit for the above Hamiltonian 
proved to have a rather interesting feature [MllTTj. At a finite time Tc = \n{gi/g2)/ (4:6), 
a pure energy crossing between two of the adiabatic states is observed. This occurs in 
the vicinity of a temporal degeneracy and results from the absence of coupling between 
the two degenerate states. As a consequence of this effect, the evolution matrix has a 
simple form, leading to a conditional entanglement between the atoms and the cavity 
mode. Furthermore, in the case of equal interactions, gi = g2, ^ number of applications 
in the field of Quantum Computing such as logic gates, state mapping and teleportation. 
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can be realised where the output can be tuned by means of a single dynamical phase 
parameter which defines the entire system evolution. 

2.1. Frequency chirps 

An interesting problem that arises when considering applications based on this system, 
is that of the fine-tuning of the output state. In this paper, we show how the Stark- 
shift technique, which is used to tune the interaction time in experiments with Rydberg 
atoms [5j, can be used to deliver a robust control over the output state. 

The whole idea requires the use of smooth Gaussian frequency chirps, such that 
each atom experiences a detuning Aj(r) at different times: 

Ai(r) = Aoexp(-(r + ro)Va2), 

(4) 

A2(r) = -Aoexp(-(r-ro)V^,2). 

These atomic detunings can be produced with Gaussian EM pulses. When the first atom 
enters the cavity, a weak EM pulse with a spatial distribution Ef{x), figured], is used to 
Stark-shift the atom, while crossing through the narrow region of Ef{x). Once the first 
atom has crossed that region, the EM pulse is turned off, long before the second atom 
enters the cavity. In a similar way, a second pulse £"1^ (x) is used after the first atom 
exits the cavity to shift the atomic transition frequency of the second atom, figure [TJ If 
both pulses have a width 2L and they have a peak at a; = —Xq and x = Xq respectively, 
then the time tq and the width as will be 
v^t + xo L 
2wo Wo 

In the adiabatic limit, the system evolution will be described by the time-dependent 
eigenfunctions of the Hamiltonian ([T]). Although the analytic expressions for the 
adiabatic states can be derived, the whole process is lengthy [16] and beyond the scope 
of this paper. Here we will demonstrate the basic features of the system by means of 
numerical simulations with the Schrodinger equation and along with simple qualitative 
arguments explain the main features of the system. 

Assuming that the Stark-chirps have a short duration, ^ 1 (L ^ wq), and 
that the location of them, r = ±ro, is away from the crossing point for the resonant 
limit, the propagator for symmetric interactions, gi ^ g2 = go, will have a similar form 
as in the resonant limit |T6l, i.e. 



|n; 61,62) |ra;6i,62), (6a) 

\n + l;gi,e2) ^ -\n + l;ei,g2), (66) 

\n + 1; 6i,5f2) cos (j)n^ \n + 1; 5-1, 62) - isin (j)n^ \n + 2; 5(1, 5(2), (6c) 

\n + 2; 5(1, 5(2) -isin \n + l;gi,e2) + cos |n + 2; 5(1, 5(2), (6d) 
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where in terms of the resonant dynamical phase the corresponding off- resonant 
adiabatic phase 0„ reads 

0„ = 0„ + 2(t/" dr(y/Ai(r)2 + A{n + 2)r/2(r) - 2T^^{T)V^r+2\ . (7) 

J — oo 

The dynamical phase is the integral over all time for the time-dependent energy of 
one of the adiabatic states in the resonant limit [T6] . 

The second term in Equation ([7]) is derived after taking into account the fact that 
when Ai(r) or A2(r) is on, then 772(1") = or rjiij) = respectively. Thus the system 
will correspond to a single atom interacting with the cavity mode, with a time-dependent 
detuning between the atomic transition and the mode frequency. Because of the choice 
made for the chirp, the adiabatic energies at this stage of the evolution will be pushed 
away from each other, and the system will adiabatically evolve during the time interval 
for which the atoms go through the fields Ej{x). 

When the chirp is off, the atom has crossed the region of the field Ej{x), the system 
returns to the initial two-atom resonant adiabatic state, with a phase shift emerging 
from the chirped interaction of each atom with the mode. Since Ai(r) = — A2(— t) and 
r7i(r) = ri2{—T) the two contributions, one for each atom, can be merged into one term 
as in Equation (JTj). Between the two chirps the system evolves as in the resonant limit, 
resulting in an additional phase first term in Equation ([7]). 

Thus in the adiabatic limit and for symmetric interactions, the system evolves 
according to ( l6al) -(l6dl). For this to be the case the previously derived conditions for 
using the adiabatic approximation must be satisfied [16]. More specifically the coupling 
between the atoms and the cavity mode must be strong, go ^ v/wq, and the delay time 
between the atoms must be of the order of Awq/v. In addition to these conditions , the 
EM fields used to produce the Stark-chirps must be far from the cavity center, xq S> wo, 
and their width must be smaller than the mode width, L <^ wq. It is also important to 
note that the pulses should be weaker than the atom-cavity coupling i.e. Aq <^ go- As 
we see later on, violation of this condition is not detrimental for the system. 

2.2. Maximally entangled atoms 

Equations fl6aj) to fl6dl) describe the conditional entanglement of the second atom with 
the cavity mode. In order for this to be the case, the atom must be initially in its ground 
state. Instead of entangling the cavity mode with one of the atoms, one could generate 
a maximally entangled state of the two atoms. Preparing the system in the factored 
state 

lV') = ^|0)(k)+ei))(b) + |e2)) (8) 

before sending the atoms through the cavity, will result in a maximally entangled state 
of the atoms 

l^en) = l\92) i\9i) - \ei)) + ^|e2) + |ei)) , (9) 
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Figure 2. The fidelity for a maximally entangled state of the atoms (O as a function of 
the chirp amplitude Aq. The fidelity was calculated after numerical integration of the 
Schrodinger equation. The parameters are: go = 30/cr, 5 = 1.25, tq = 2.0 and = 0.2. 
The first peak where a maximally entangled state is formed is for Aq ~ QAAqq. The 
value obtained from Equations ([7]) and ((TO)) is Aq ~ QAQg^. 



if (pn = '2,17171 with m an integer. Thus, a maximally entangled state of the two atoms 
can be generated and the output can be fine-tuned by means of the Stark-chirps. 

In figure [2] the fidelity for a maximally entangled state [3] is given for different values 
of the chirp amplitude Aq. The first important feature is the periodic reappearance of a 
maximally entangled state even for Aq > Qq. This is due to the fact that the system still 
evolves adiabatically, even if the Stark fields are stronger than the atom-cavity coupling. 
For Aq < go the phase 0„ oc Aq and the corresponding error for the fidelity 



will be a linear function of Aq. Thus the Stark-shift technique is expected to be robust 
for weak chirps Aq ^ (^o- 

The robustness is demonstrated in figure [3] where the variations for the fidelity 
around the optimum value A^ax = 0.44(7o are plotted as a function of the ratio Ao/A^ax- 
From this we can see that variations of the order of 10% around the optimum value 
Aq = ^max havc a very small impact on the fidelity, (< 1%). Thus generating a 
maximally entangled state, or reahsing the previously proposed applications [16l[I7j, can 
be robust by means of weak EM fields used to detune the atomic transition frequencies 
from the cavity frequency. This in return will induce a second order shift in the 
dynamical phases ([7]), which can be used to fine-tune the output of the evolution under 
the Hamiltonian ([1]) in the adiabatic limit. 




(10) 
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Figure 3. The fidelity for a maximally entangled state (solid) for variations of Aq 
around Amaa; = 0.44(70- The fidelity (dashed) for a high quality cavity, Q = 4.2 x 10^° 
|20j . and circular Rydberg atoms with lifetimes Tat = 30ms (5j. Other parameters: 
50 = 30/ct, 5 = 1.25, To = 2.0 and ct, = 0.2. 

3. Entanglement spatial properties 

Up to this point, one of tlie main assumptions was that both atoms are coupled to the 
cavity field via time- dependent coupling functions with equal amplitudes gi = g2- In 
the most general case and for a spherical mirror resonator [TT] , these amplitudes will be 
functions of the coordinates yj and Zj, figured! of the form 

djiVj^Zj) = gocos{kzj)exp (-?/|/(2wo)^) , (11) 

where A; = 27r/A is the mode wavenumber. From this we see that gi ~ g2 if Zi ~ Z2 and 
Vi ~ 1/2- What is now interesting, is the properties of entanglement when this latter 
condition is not satisfied. 

For the resonant limit Aj = 0, asymmetries in the coupling profiles ([2]) would 
introduce a second adiabatic phase 9n, [E], in addition to (pn- Despite this the 
conditional entanglement between the second atom and the mode still remains, but 
now the output state could have two different forms each defined by one of the two 
angles 6'„ or This qualitative picture holds even when the chirps (jlj) are taken into 
account. An important feature, is that 6'„ unlike 0„, does not experience any shift as 
0„ does, see Equation (JTl). 

It is now interesting to consider the entanglement properties for the two atoms, as 
their positions Zj are varied. In order to take into account the possibility of forming a 
tripartite entangled state of the atoms with the cavity mode, we quantify entanglement 
in terms of the multipartite pure state concurrence [23], which for our case is 

Cs im = ^3-trpf-trpi-trp2, (12) 

where pj is the reduced density matrix for atom j and pc is the reduced density matrix 
for the cavity mode. Furthermore, we assume that the first atom moves along an anti- 
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Figure 4. The concurrence as a function of the atomic separation Z2 — zi. The 
initial state is 1-0) The parameters are: go = SO/u, 5 — 1.25, Aq — 0.44go, tq — 2 
and (Ts = 0.2. 



node of the standing wave, i.e. zi = mA, while for the second atom Z2 is arbitrary. With 
these conditions the couphng strengths will be 

9i = 9o, g2 = go cos {k{z2 - zi)) . (13) 

In figure m we plot the results of a numerical simulation for the concurrence ( IT2ll . 
as a function of the atomic separation 2:2 — ^1. As one expects for integer or half integer 
values of the ratio {z2 — zi)/X we have that gi = ±(72, and a maximally entangled state 
of the two atoms is formed. These points are the local minima where C3 (lip)) = 1. 
On the other hand, whenever the atomic separation is {z2 — Zi)/X = (2m + l)/4, then 
g2 = and this corresponds to the single atom Jaynes-Cummings model. For this limit, 
the first atom entangles with the cavity mode, forming sharp dips in figure H] where the 
concurrence obtains its minimum value. Around the local minima where a maximally 
entangled state of the atoms is formed, symmetric spikes with C3 (lip)) > 1 signify the 
tripartite entanglement of the atoms and the cavity mode. For values of C3 (lip)) < 1 
and arbitrary atomic separations, the system could be in either a tripartite or bipartite 
entangled state. In order to distinguish between the two, additional information is 
necessary such as the values of tr pj and tr pi or the populations and phases for the bare 
states. 



4. Spontaneous emission and cavity losses 



Spontaneous emission from the atoms to modes other than the cavity mode, and photon 
losses through the cavity mirrors are both detrimental for entanglement. In an attempt 
to quantify the importance of both effects, we solved the master equation for the density 
matrix of the entire system, 
dp 



dt 



— 1 



H{t),p +Cs{p) + C,{p). 



(14) 
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The Liouvillian terms Cs{p) and Cc{p), respectively describe the atomic spontaneous 
emission with a rate F, and the decay of the cavity with a rate 7 into different thermal 
reservoirs at zero temperature |24j . 



i=i-2 (15) 
^c{p) = —— (a^ap + pa) a — 2apa)^ . 

Calculating multipartite concurrence for a mixed state is not an easy task [25] . One 
would have to perform an optimization for the concurrence C3 (iV'i)) (!T2l) over all pure 
state ensembles which equally represent the mixed state p that corresponds to 

the solution of Equation ( |T^ . Instead of this approach, we were able to show after a 
number of simulations, that the final state of the system has a simple form 

p(oo) ^ |0)(0|®p,(oo), (16) 

where Pa(oo) is the reduced density matrix for the atoms for t = 00. This result is for 
the initial state ([8]), and is rather accurate since the average photon number {n{oo)) is 
very small, (n(oo)) < 10^'^, whereas the correlations between the mode and the atoms 
are negligible. This result, allows the calculation of the concurrence for the atomic pair, 
by means of the two-qubit mixed state concurrence [261 EZ] 

c(p) = max |o, V^i — "\/A2 — V^s — V^4| , (17) 

where Xj are the eigenvalues of the matrix 

R{p) = p{(Xy®ay)p*{ay^ay), (18) 

where ay is the Pauli matrix [28j and the eigenvalues Aj are in increasing order i.e. 
Ai > A2 > A3 > A4. 

The results from numerical simulations with Equation f[T^ and (1171) are plotted in 
figure O where we see that entanglement between the two atoms exponentially decays 
with respect to both decay rates. Furthermore, the most important detrimental effect 
is that of spontaneous emission from the atoms. This is because the effective decay rate 
for the subsystem of the two atoms is greater than the single atom decay rate F. That 
is, entanglement is more sensitive to decay due to spontaneous emission than decay due 
to cavity losses. 



4.1. Experimental feasibility 

For experiments with circular Rydberg states of Rb^^ atoms, with lifetimes Tat = 30ms 
[5l [21] , the main requirement for adiabatic evolution is for the system to be in the strong 
coupling regime, QoWq/v ~ 20 — 30. For a coupling strength go/^-K = 50kHz and a mode 
waist 2wo = 6mm, the atomic velocities must be f ~ 60 — 95m/s. These speeds are 
large enough to avoid reflection of the atoms from the cavity field. In addition to this 
the initial atomic displacement must be 2vAt (4 — 5)wq = (12 — 15)mm. 



Entanglement in adiahatic cavity QED 



10 




I ' ' — r 

0.05 0.1 0.15 0.2 



r(l/a) 

Figure 5. The two-qubit concurrence P7)l . for different cavity decay rates 7 with 
r = (top) and for different spontaneous emission rates T with 7 = (bottom). 
Other parameters: go « 18.9286/(7, 5 = 1.25 and Aq = 0. The initial state is ([8]), 
whereas for this parameters and for 7 = F — the output is a maximally entangled 
state ©. 



This distance is bigger than the width of the Stark fields which is 2L ^ OAwq = 
1.2mm. This allows the Stark fields to be turned on and off at proper times so that 
they interact with only one atom and consequently produce the desired chirps dl]). 
The location of the Stark pulses can vary between 4.5 and 6mm with no impact on 
the robustness of the system. In addition to all these, a recent development was the 
engineering of a cavity with Q = 4.2 x 10^*^ [20l [2T] . For this quality factor and the 
corresponding lifetimes for the circular Rydberg atoms, the fidelity for a maximally 
entangled state ([9]) is reduced by about 3% varying between 96 and 97%, figure [3], which 
is a very high fidelity. 

5. Conclusion 

In this work we have studied a time-dependent, off- resonant two-atom Tavis-Cummings 
Hamiltonian. Considering pairs of slowly moving atoms, sequentially crossing a 
Gaussian mode resonator, we examine the adiabatic limit for the system. As we were 
able to show, for proper choice of atomic frequency chirps the evolution of the system 
bears strong similarities to the resonant limit. The system propagator has the same 
structure as in the resonant limit, and one can realise the exact same applications 

[mini. 

The only and also crucial difference from the resonant limit, is that the adiabatic 
phase is now shifted due to frequency chirps. Because of the linear dependence of this 
shift with respect to the chirps amplitude, all the proposed applications can be robust. 
As an example, the creation of a maximally entangled state of the atoms can be achieved 
with fidelities higher than 99%. Despite potential detrimental effects, such as atomic 
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spontaneous emission or photon losses from the resonator mirrors, recent developments 
in the fabrication of high quality cavities [201 121], allow fidelities as high as 97%. In 
addition to this, the adiabatic condition can be easily satisfied for experiments with 
Rydberg atoms in the microwave regime [5]. 

Although time- dependent frequency chirps could be detrimental for the adiabatic 
approximation, our results show that proper timing of these chirps could enable 
adiabatic evolution even if they are stronger than the coupling between the atoms and 
the field. One could again coherently tune the output of the evolution, but this would 
be less robust. This is due to the sinusoidal dependence of the error with respect to the 
chirp amplitude for strong Stark-fields. On the other hand, for weak chirps the error in 
the fidelity is very small making the Stark-shift technique highly robust. 
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